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Abstract. Single diffractive Drell-Yan reaction in hadron-hadron collisions is considered as 
an important source of information on the properties of soft QCD interactions. In particular, 
it provides an access to the dynamics of the QCD factorisation breaking due to the interplay 
between hard and soft interactions which leads to a nontrivial energy and scale dependence of 
the Drell-Yan observables. We study the process at forward rapidities in high energy proton- 
(anti)proton collisions in the color dipole approach. Predictions for the total and differential 
cross sections of the diffractive lepton pair production are given at different energies. 



1. Introduction 

Diffractive large rapidity gap processes in QCD constitute a noticeable fraction of all observed 
events. Depending on energy, it may vary from a fraction of percent in pp up to 10 % and more 
in ep scattering. Generally, such processes are very sensitive to the soft and nonperturbative 
interactions despite the presence of a hard scale, which makes them extremely difficult to 
investigate from both the theoretical QCD and experimental points of view, and intrinsic 
uncertainties of existing phenomenological models are still quite large. 

Diffractive Drell-Yan process at forward rapidities [T], as well as diffractive heavy flavor 
production [2], is one of such processes, which give us an immediate access to soft QCD evolution 
close to saturation regime. The understanding of the mechanisms of inelastic diffraction came 
with the pioneering works of Glauber [3], Feinberg and Pomeranchuk [3], Good and Walker 
[5]. If the incoming plane wave contains components interacting differently with the target, the 
outgoing wave will have a different composition, i.e. besides elastic scattering a new diffractive 
state will be created (for a detailed review on QCD diffraction, see Ref. [6]). In our case, such 
a new state is given by the deeply virtual photon radiation in the forward direction, which then 
can be seen as e.g. heavy fJ,~^iJ.~ pair in a detector. 

The single-diffractive Drell-Yan reaction in pp collisions is characterized by a relatively small 
momentum transfer between the colliding protons, such that one of them, e.g. pi, radiates a 



hard virtual photon k'^ = ^ and hadronizes into a hadronic system X both moving 
in forward direction and separated by a large rapidity gap from the second proton p2, which 
remains intact, i.e. 

Pi+P2^-f*{l^l~) + X + {gap)+p2 (1) 

Both the di-lepton and X stay in the forward fragmentation region. In this case, the virtual 
photon is predominantly emitted by the valence quarks of the proton pi. We will refer to this 
as the diffr active Drell-Yan process at forward rapidities. 

The dipole approach, previously applied to diffractive Drell-Yan reaction in Ref. [1], led to the 
QCD factorisation breaking, which manifests itself in specific features like a significant damping 
of the cross section at high ^/s compared to the inclusive DY case. This is rather unusual, since 
a diffractive cross section, which is proportional to the dipole cross section squared, could be 
expected to rise with energy steeper than the total inclusive cross section, like it occurs in the 
diffractive DIS process. At the same time, the ratio of the DDY to DY cross sections was found 
in Ref. [Ij to rise with the hard scale, M^. This is also in variance with diffraction in DIS, which 
is associated with the soft interactions [7J. 

Long-range soft interactions between target and projectile particles, or the absorptive 
corrections affect differently the diagonal and off-diagonal terms in the hadronic current |8j, in 
opposite directions, leading to an unavoidable breakdown of the QCD factorisation in processes 
with off-diagonal contributions only. Namely, the absorptive corrections enhance the diagonal 
terms at larger ^/s, whereas they strongly suppress the off-diagonal ones. In the diffractive 
DY process a new state, the heavy lepton pair, is produced, hence, the whole process is of 
entirely off-diagonal nature, whereas in the diffractive DIS contains both diagonal and off- 
diagonal contributions [6]. 

Another reason of the QCD factorisation breaking is more specific and concerns the interplay 
of soft and hard interactions in the DDY amplitude. In particular, this leads to the leading 
twist nature of the DDY process, whereas DDIS is of the higher twist [IJ. Large and small 
size projectile fluctuations contribute to the diffractive DY process at the same footing, which 
further deepens the dramatic breakdown of the QCD factorisation in DDY. This work is devoted 
to a detailed study of consequences of such a breakdown in typical DY observables. 

2. Diffractive Drell-Yan in dipole-target scattering 

The hard part of the Drell-Yan process is given by the inelastic amplitude of 7* radiation by a 
projectile quark (valence or sea) due to its interaction with the target through a gluon exchange 
as shown in Fig. [TJ It consists of two terms corresponding to interaction of two different Fock 
states with the target - a bare quark before the photon emission \q) (s-channel diagram), and a 
quark accompanied by a Weizacker- Williams photon |(?7*) (tt-channel diagram). 

Elastic dipole scattering depicted in Fig. [1] corresponds to forward scattering at small 
momentum transfers in the t-channe0. In the leading order, the elastic scattering amplitude 
is given by one- loop diagram with two t-channel gluon exchanges. For on-shell intermediate 
spectators, corresponding four-dimensional loop integral can be reduced to two-dimensional one 
over the transverse momentum of one of the gluons 

2iImF,i{6^) = J ^ A{k^)A{6^ - k^), 6^ « \k^\ (2) 

where A represents (inelastic) amplitude for one t-channel gluon exchange, and the last strong 
inequality guarantees that the proton target survives the scattering, hence, the elastic nature of 

^ Generally speaking, \/—t — >■ Aqcd corresponds to the physical forward scattering limit since transverse 
momentum of a proton in the final state cannot be resolved to a better accuracy than its inverse size. 
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Figure 1. Leading order contribution to the diffr active Drell-Yan in the dipole-target colHsion. 



the process. Then the convolution theorem of Fourier analysis leads to the optical theorem 



2ilmU{b) = \A{b)\' 



(3) 



In the case of multiple elastic rescattering, this relation leads to eikonalization of the bare elastic 
amplitude which then correctly takes into account gap survival effects. However, if one uses the 
elastic amplitude f^i fitted to soft data, it must already contain the soft unitarity corrections, 
and one should include them twice as is sometimes done in the literature. 

Straightforward calculations lead to qq dipole scattering amplitudes for s and n-channel 
photon emission, respectively. 
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where the last terms subtract the contributions from diagrams corresponding to the situation 
when none of the gluons couple to the same quark line with the hard photon. Then, implied the 
fact that all fields disappear at infinite separations, i.e. Vq{b) — >• when |6| — )■ oo, we have due 
to antisymmetry of the integrand 



yg(6 + r;)|2"- |yg(6 + r; + af)|^' 
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(4) 



such that these terms do not contribute to the final result. Using the optical theorem for the 
elastic amplitude 



2ilm Ub, fp) = j^J2Il \^i(b) - Vgib + 

X CfCi 



we can finally write 



Mkr,' [0, rp, r, a) 



2Im fei{b, fp) - 2Im fei{b, fp + af) 



(5) 



i.e. the amplitude of the diffractive radiation is proportional to the difference between elastic 
amplitudes of the two Fock components, with and without the photon radiation. When a quark 
fluctuates into the upper Fock quark-photon state with the transverse separation r, the final 
quark gets a transverse shift Ar = ar. Then the quark dipoles with different sizes in the |2q) 
and 12(77*) components interact differently, and their difference corresponds to the diffractive 
Drell-Yan amplitude ([5]). 



3. Diffractive Drell-Yan in proton-proton scattering 

In the dipole picture, the typical enhanced Regge graphs correspond to elastic scattering of 
higher Fock states, which contain gluons, e.g. \qqg), \qqgg), etc. Note that in our approach 
we take into account the lowest Fock state \qq) contribution only. Such an approximation is 
justified for not very small fraction Xyi = p'i^/pf and scale M^, where valence/sea quarks are 
dominated and the gluon contribution is rather small. 

The total hadronic amplitude of the diffractive Drell-Yan process can be written as [1] 

A.f=A^+A^+A!^, (6) 

where each term corresponds to 7* radiation by one of the valence (sea) quarks in the proton, 
in particular, 

A^^j) {x^i,q±, X^) = ^ Q!^ y Cpri(fr2(f'r^d?r(fihdxq^dxq2dxq^ 

X '^i{ri,r2,r^]Xq^,Xq2,Xq,^)^*f{ri + ar,r2,r:i-Xq^ - x^i,Xq2,Xq^) 

Mg-; ^""i - ^""2' ^ «) + n - ^""s, r, a)] e^'^^-V^^^-? (7) 

Here, = L,T; l± = 5^ — q^ja (z-axis is directed along initial proton momentum); the hard 
photon with virtuality q^ = ^ m^, transverse gj_ and fractional longitudinal x^i momenta 
is emitted from the first valence quark with impact parameter ri, other two valence quarks in 
the proton have impact parameters r2 and rs, respectively; r is transverse separation between 
the photon and the radiating quark; a = x^i jxq^ is the fraction of longitudinal momenta taken 
away by the photon from the radiating quark; Mqq are the Fourier-transformed amplitudes for 
the elastic quark dipole scattering off the proton target accompanied by the hard L, T-polarized 
photon emission; ^'j j are the light-cone wave functions of the 3g systems in the initial and final 
state, respectively. In Eq. ([7]) we implicitly assumed that exchanges i-channel gluons all together 
take a negligibly small longitudinal momentum compared to the collisions energy -y/i and, hence, 
corrections to quark momenta due to gluon couplings are neglected in the wave functions. 

The differential cross section for the single diffractive di-lepton production in the target rest 
frame reads 

(n)i ^ . da d'^q±d'^5± 



dM^dn 
167r2M4' 



X aeme^iX^)el{X^)L'^^—TJ7I^ X^ = L,T (8) 



where prefactors provide averaging over colors and helicities of exchanged t-channel gluons, p^^ 
is the energy of the rci(iia,ting Tzth. qucirk in the initicil sta^te, ti — 1, 85 Q^em — 

eV(47r) = 1/137 

is the electromagnetic coupling constant. The second line in Eq.dH]) describes decay of 7* into the 
leptonic pair into solid angle dCl = dcpd cos 9, and L^^'^ is the standard leptonic tensor. We 
keep in the cross section only diagonal in the photon polarization Xy = L,T terms (non-diagonal 



ones drop out after integration over leptonic azimuthal angle 0). Integrating the diffractive 
differential DY cross section over the photon transverse momentum q±_ we get 



dlnadW^ <P5^ 
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we get the diffractive 7* production cross section in the following differential form 
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2Im fi - f2) + 2Im ff^^{h, f 1 - f 2 + af) 
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where the Kopeliovich-Schafer-Tarasov (KST) parameterization of the elastic dipole-target 



amplitude fitted to the soft data O [TOl [TT] and, hence, valid at |rj 
Finally, going over to the forward limit = we obtain 



b, i ^ j is used. 
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We see that normalization of the cross section agrees with the original result of Ref. [T]. The 
total diffractive cross section is then given by 



da{pp P"y*X) 
din a 



1 d^a{pp — > p'y*X) 
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(14) 



where [s) is the diffractive slope similar to the one measured in diffractive DIS. 

The next step is to introduce the proton wave function assuming the Gaussian shape for the 
quark distributions in the proton as 
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where a = (r|g) and b = {Rq) are the diquark mean radius squared and the quark mean distance 
from the diquark squared, respectively. In this work, we will use the simplest case of symmetric 
valence quarks distribution assuming that Vqq = Rq = 0.85 fm. 



Then valence quark distribution in the proton is given by 

dXq^dXq^pi^Xq^, Xq2, Xq^) = Pqi{Xqi^ 



where we integrated out the longitudinal fractions of the diquark in the proton. Generalization 
of the three-body proton wave function (jlSp including different quark and antiquark flavors, as 
well as sea quarks, leads to the proton structure function [12] 



^qiPgi^l) + Pqi^q)] = ^I^2{Xq), Xq = ^ . (16) 

q Xq a 

4. Numerical results for differential DDY cross sections 

In Fig. [2] the ratio of the diffractive to inclusive DY cross sections is plotted as a function of 
di-lepton invariant mass squared (left panel) and photon fractional light-cone momentum 
x^i (right panel) at different energies. In the left panel, the curves are given for fixed x^i = 0.5 
(solid lines) and x^i = 0.9 (dashed lines). In the right panel, the curves are given for fixed 
M2 = SOGeV^ (sohd lines) and = 500 GeV^ for ^ = 14 TeV and 500 GeV, and 
M2 = 200 GeV^ at = 40 GeV (dashed hues). The pairs of solid/dashed curves in the 
both panels correspond to y/s = 40 GeV, 500 GeV and 14 TeV from top to bottom, respectively. 
Here we used the KST parameterization for the dipole-target scattering amplitude [10\ [9l [11] 
and F2 parameterization by Cudell and Soyez [13] are used here. In this calculation we consider 
the unpolarized case summing up the contributions of longitudinal and transverse parts both in 
the diffractive and inclusive cross sections. 
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Figure 2. The ratio of the diffractive to inclusive Drell-Yan cross sections as function of the lepton-pair 
invariant mass squared (left panel) and photon fraction x-yi (right panel) at different energies. 



As seen from Fig. [21 the DDY-to-DY cross section ratio is falling with energy. However, 
naively one could expect basing on QCD factorisation, that the DDY cross section, which is 
proportional to the dipole cross section squared, should rise with energy steeper than the total 
inclusive cross section. At the same time, the ratio rises with the hard scale of the process, M^. 
This also looks counterintuitive, since diffraction is usually associated with soft interactions |14| . 
These effects are different from ones emerging in Regge factorisation-based calculations, where 
we observe a slow rise of the DDY-to-DY cross section ratio with c.m.s. energy and its practical 
independence on the hard scale of the process |15J. 

In order to understand such an interesting behavior of the DDY-to-DY cross sections ratio 
obtained in the color dipole approach, let us look at the amplitude of the DDY process, which 



is proportional to the difference between the dipole cross sections of the Fock states with and 
without the hard photon emission [1], i.e. 

Mflw~fT(i?)-a(i?-af) = ^^e-f^'/-f^o(^) (f-R)+h.o. (17) 

assuming the simplest Golec-Biernat-Wusthoff (GBW) slope for the dipole cross section [Mj, and 
the hardness of the emitted photon implies r ~ 1/M ^ Rq{x). We see now that the diffractive 
DY amplitude is linear in r, so the diffractive cross section turns out to be a quadratic function of 
r, which is different from e.g. the diffractive DIS process where the cross section is proportional 
to and is dominated by soft fluctuations (see e.g. Refs. [HIT]). Since the diffractive DY cross 
section is proportional to r^, then soft and hard interactions contribute on the same footing [i], 
which is one of the basic sources of the QCD factorisation breaking in diffractive DY process. 




Figure 3. Diffractive Drell-Yan cross section (in fb) as function of the lepton-pair invariant mass 
squared (left panel) and photon fraction x^i (right panel) at the RHIC II c.m.s. energy ^/s = 500 
GeV. 



We also compare predictions for the diffractive DY cross section for different 
parameterizations for elastic dipole-target scattering amplitude corresponding to scattering of 
small (GBW given by Refs. [ISlEllIZj) and large (KST given by Refs. 01101 HJ) dipoles. As an 
example, in Fig. [3] we present the diffractive Drell-Yan cross section as function of the lepton- 
pair invariant mass squared M'^ (left panel) and photon fraction Xji (right panel) at the RHIC 
II c.m.s. energy ^/s = 500 GeV. We notice that the GBW parameterization leads to roughly 
a factor of two smaller cross section than the one obtained with the KST parameterization, 
however, both of them exhibit basically the same x-yi and shapes. It means that the 
evolution of the dipole size can only affect the overall normalization of the DDY cross section. 
Since arguments in the elastic amplitude fei{rp,b), the impact distance between the target and 
the projectile b and the transverse distance between projectile quarks rp ^ \ fi — fj\, i ^ j, are 
of the same order and given at the soft hadronic scale, then the use of KST parameterization 
fitted to the soft hadron scattering data data is justified in the case of diffractive DY. 



5. Conclusion 

The QCD factorisation breaking effects in the diffractive Drell-Yan process lead to quite 
different properties of the corresponding observables with respect to QCD factorisation-based 
calculations. A quark cannot diffractively radiate a photon in the forward direction, whereas 
a hadron can due to the presence of transverse motion of spectator quarks in the projectile 
hadron. For this reason, the diffractive DY cross section depends on the hadronic size explicitly 
breaking the QCD factorisation. 



This leads to the physical picture where hard and soft interactions are equally important 
for DY diffraction, and their relative contributions are independent of the hard scale, like in 
the inclusive DY process. This is a result of the specific property of DY diffraction: its cross 
section is a linear, rather than quadratic function of the dipole cross section. On the contrary, 
diffractive DIS is predominantly a soft process, because its cross section is proportional to the 
dipole cross section squared. 

Contrary to what follows from the calculations based on QCD factorisation, the ratio of the 
diffractive to inclusive cross sections falls with energy, but rises with the di-lepton effective mass 
M. This happens due to the saturated behavior of the dipole cross section which levels off at 
large separations. All these properties are different from those in the diffractive DIS process, 
where QCD factorisation is exact. In addition, we made predictions for the differential (in 
photon fractional momentum x-yi and di-lepton invariant mass squared Af^) cross sections for 
the diffractive DY process at the energies of RHIC (500 GeV) and LHC (14 TeV). 
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